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INTRODUCTION 


High  resolution  sonar  systems  are  used  for  bottom  mapping,  small 
object  location,  and  other  underwater  tasks  requiring  high  quality 
images.  The  term  high  resolution  sonar  is  used  here  to  describe  sys- 
tems operating  in  the  frequency  range  of  100  kHz  and  greater,  with  beam- 
widths  of  less  than  2 degrees  and  range  resolution  of  better  than  1 
metre.  Unfortunately,  most  theoretical,  high  resolution  sonar  design 
and  performance  evaluations  are  based  on  solution  of  the  "sonar  equa- 
tion," and  many  important  phenomena  are  not  considered  when  only  energy 
relationships  are  evaluated.  Examples  are:  the  relative  importance  of 

sidelobe  levels  and  sonar  beamwidth,  display  dynamic  range,  and  range 
resolution. 

The  work  presented  here  is  the  result  of  a preliminary  investiga- 
tion of  advanced  simulation  of  high  resolution  sonar  systems  and  the 
operating  environment.  The  purpose  of  this  basic  model  is  to  demon- 
strate the  feasibility  and  versatility  of  such  a computer  simulation 
and  to  provide  a foundation  for  a more  detailed  model. 

The  major  elements  of  the  model  are  the  target  area,  the  projector- 
hydrophone  arrays,  the  beamformer,  and  the  display.  A schematic  of  the 
simulation  concept  is  shown  in  Figure  1.  The  target  area  is  modeled  by 
defining  a reflectivity  amplitude  and  phase  for  incremental  elements 
within  the  target  region.  The  projector  array  is  modeled  by  describing 
a line  array  length  in  coordinates  of  the  array  center.  The  hydrophone 
array  is  modeled  by  a number  of  discrete  elements  of  finite  length  and 
a three  dimensional  coordinate  for  each  element.  The  beamformer  is 
modeled  as  either  an  FFT  algorithm  using  outputs  of  the  discrete  hydro- 
phone elements  as  inputs,  or  as  a convolution  of  the  hydrophone  array 
spatial  impulse  response  with  a target  reflectivity  distribution.  The 
display  is  driven  by  a mini-computer  which  can  be  programmed  to  scale, 
expand,  or  compress  image  data  as  desired. 


TARGET  ECHO  SIMULATION 


The  geometry  of  the  target  and  array  simulation  is  shown  in  Figure  2. 
The  target  area  is  modeled  by  defining  a reflectivity  amplitude  and  phase 

(Text  Continued  on  Page  3) 
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FIGURE  1.  BASIC  ELEMENTS  OF  SONAR  SIMULATION 
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FIGURE  2.  GEOMETRY  OF  SONAR  TARGET  AND  ARRAY  SIMULATION 
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for  incremental  elements  within  the  target  region.  These  arrays  are 
designated  TA(m,n)  and  Tp(m,n)  for  the  amplitude  and  phase  components, 
where  (m,n)  represent  coordinate  points  in  the  target  area.  Empirical 
data  were  used  to  establish  echo  levels  from  various  bottom  types  as  a 
function  of  frequency,  grazing  angle,  and  bottom  type  such  as  sand,  mud, 
etc.;  however,  the  statistics  of  this  reverberation  are  not  well  defined. 
Therefore,  a random  number  procedure  with  only  upper  and  lower  bounds 
specified  was  used  to  fill  the  target  arrays  with  the  background  ampli- 
tude. Background  amplitude  numbers  were  in  the  range  0.05  to  0.2  and  the 
target  amplitude  numbers  in  the  range  0.5  to  1.0.  Phases  were  taken  as 
constant  or  random  from  0 to  i. 


The  projector  is  modeled  by  assuming  a line  array  length  L posi- 
tioned at  a point  (0,0,h).  The  directional  response  of  the  projector 
array  is  computed  by 


PA<0) 


sin(ir(L  /X)  sin  0)/(ir(L  /X)  sin  0) 
P P 


P (0)  = sgn(P. (0) ) (i: 

P A 

where  0 is  measured  from  the  origin  with  respect  to  the  y axis,  and  X 
is  the  wavelength.  The  letters  sgn  (P^(0)  denote  the  signum  function 
and  is  +1  if  PA(0)  1.  -1  if  Pa(0)  < using  Equation  (1),  the 

assumption  is  made  that  the  target  area  is  in  the  projector  far  field; 

»max  > Lp  /X,  where  range  R is  measured  from  the  origin  and  is 
given  by  R = (x2  + y2)^.  The  complex  amplitude,  At,  reflected  from  a 
target  point  as  measured  in  the  plane  of  the  target,  is  given  by 


A (m,n)  = 


Vm’n]  I PA(0) 


where  0 is  the  slant  range  from  projector  center  to  target  point,  a is 
absorption  coefficient  in  dB/unit  distance,  and  S is  the  projector  source 
level  in  dB. 


TARGET  ECHO  PROCESSING 


Echo  processing  simulation  includes  sonar  beamforming  and  the  effects 
of  a finite  pulse  length.  Included  are  the  parameters  of  the  hydrophone, 
aperture  apodlzation  or  shading,  FFT  characteristics,  pulse  shape,  and 
display  characteristics. 

The  receive  hydrophone  is  modeled  by  assuming  k number  of  elements 
each  of  receiving  length  and  center-to-center  separation,  s^.  The 
total  length  of  the  hydrophone  is  given  by  = k ■ s^  + t,^. 
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Discrete  elements  of  finite  length  are  assumed  so  the  reflected 
wavefront  may  be  properly  sampled  and  beamforming  may  be  simulated  using 
a Fourier  Transform  algorithm.  Focusing  may  be  accomplished  by  placement 
of  element  coordinates  to  achieve  array  curvature. 

The  complex  amplitude,  (n,j),  output  of  a discrete  element  i re- 
sulting from  returns  from  reflectors  in  target  row  n may  be  computed  as 


mt 

= £ Attm>n]  10 

m=l 


-iCf1  r ) (3) 

A m,n,j 


where  rm  n>j  is  the  slant  range  from  target  element  (m,n)  to  hydrophone 
element  j . ’ 

In  many  systems,  aperture  apodization  or  shading  is  used  to  reduce 
sidelobe  (first  diffraction  order)  levels^1^.  Shading  is  accomplished 
normally  by  amplitude  shading,  but  any  complex  operation  on  the  outputs 
of  discrete  elements  may  be  simulated  by  multiplying  A^(j)  by  the  appro- 
priate complex  function  S^(j). 

Beamforming  is  accomplished  by  the  Cooley-Tukey  FFT  algorithm^ . 

The  FFT,  of  single  dimension  N,  is  performed  using  as  input  data  one 
range  increment  n at  a time.  Using  the  values  calculated  for  discrete 
transducer  element  outputs  as  the  input  array,  the  operation  performed 
by  the  discrete  FFT  operation  is 


Q [n,ql  = £ A.  [n,j]  S , [j]  e 

j=l 


(4) 


where  Q is  a complex  number  representative  of  the  beamformer  output.  In 
the  computer  model  the  complex  value  of  Q is  contained  in  the  real  and 
imaginary  arrays  Q and  Q , each  having  dimension  n x N. 

K 1 L 

Each  element  of  the  array  represents  a sonar  beam  output  in  the 
direction  given  by  the  expression 


e 


arcsin  (q  • A/N  • S^) 


(5) 


where  q is  the  qth  array  element  in  the  FFT  output  and  takes  on  values 
-N/2  to  +N/2  . 


^R.  C.  Hansen,  Microwave  Scanning  Apertures,  Vol.  1,  Academic  Press, 
New  York  (1964) . 

(a)J.  W.  Cooley  and  J.  W.  Tukey,  Mathematics  of  Computation,  Vol.  19, 
p.  297  (1965). 
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The  beamforming  or  FFT  operation  is  carried  out  separately  for  each 
target  row,  n.  Therefore,  the  arrays  QR  and  Qj  are  two-dimensional 
where  the  parameter  n is  representative  of  range  and  q is  representative 
of  azimuth.  To  simulate  a pulse  echo  system  having  a pulse  length  ex- 
ceeding the  length  of  incremental  range  elements,  it  is  necessary  to 
compute  a running  sum  or  convolution  on  Q,  with  a rectangular  function 
representing  the  length  of  the  pulse.  If  a pulse  of  length  v incremental 
range  units  is  to  be  simulated,  the  values  contained  in  the  Q arrays  are 
modified  as  follows 

n+v 

Qf  [n,q]  = £ Q [ti,  q]  (6) 

n 

where  Q represents  the  results  of  the  running  sum  over  the  pulse 
length  v. 

The  image  amplitude  of  the  target  region  can  now  be  written  as 

I[n,q]  = (QRf n, q]  + Q^n.q])^  . (7) 

This  amplitude  represents  a degraded  version  of  the  target  region  re- 
flectivity which  has  been  degraded  first  by  the  impulse  response  of  the 
array,  and  secondly  by  the  finite  pulse  length.  An  image  of  the  target 
area  may  be  obtained  by  converting  the  numerical  values  in  array  I to 
density  variations  on  photographic  film  or  by  intensity  modulation  on  a 
CRT.  The  original  and  undegraded  target  data  can  be  compared  by  trans- 
lating the  numerical  values  in  array  T to  object  amplitudes  0[m,n]  by 

-t(~-  T [m,n] ) 

0[m,n]  - TAe  p (8) 

and  using  these  values  to  intensity  modulate  a CRT, 


CONVOLUTION  MODEL 


The  model  previously  described  is  suitable  for  forward-looking,  multi- 
beam sonar  systems.  In  side- looking  systems,  a single  beam  is  formed  by 
the  sonar  array  and  is  caused  to  scan  the  sea  bottom  by  forward  motion 
of  the  sonar  platform.  For  this  geometry,  a model  based  on  convolution 
is  more  suitable. 

In  this  model,  a target  array  is  defined  in  Equation  (1).  A direc- 
tional response  of  the  projector  is  calculated  and  the  complex  reflectivity 
distribution,  A^,  is  determined  from  Equation  (2).  A receiving  array 
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I directional  response,  P^O),  is  computed  from  receiving  array  geometry 

or  data  from  measured  patterns  are  used  if  they  are  available.  The  beam- 
forming process  is  described  by  the  convolution  of  P.  [m,n]  and  At_'m,n], 
where  P^tm.n]  is  the  hydrophone  array  response  transformed  to  target 
coordinates  with  the  response  curve  maximum  centered  at  m . The  convolu- 
tion is  accomplished  through  solution  of 

IQ[mo,n]  = £ Ph[m-mo,n]  (9) 

m 

where  [m0,m]  represents  the  image  of  the  point  mQ  in  the  target  area. 

This  summation  must  be  carried  out  for  each  of  the  points  (m,n)  in  the 
target  array.  The  results  of  this  calculation  are  stored  in  arrays  Qr 
and  Qj,  and  pulse  length  effects  can  be  simulated  by  Equation  (6).  The 
image  is  then  described  by  Equation  (7). 


APPLICATION  OF  SONAR  MODELS 


A simulation  was  made  of  a sonar  target  configuration  to  demonstrate 
the  results  of  applying  the  FFT  model: 

m = 60,  nt  = 40,  N = 128 

H = 10  metres,  R . =20  metres,  R =60  metres 

min  max 

L = 20X 
P 

s^  = 0.5X,  8^  * 0.1X,  k,v  = variable 

0.5  < T^  < 0.2,  0 < Tp  < it  (background  parameters) 

T.  = 1,  T =0  (target  parameters) 

A p 

For  the  simulation  results  shown  in  Figure  3,  a target  of  8 x 8 
metres  was  simulated  by  a reflectivity  value  of  unity  with  a constant 
phase  angle.  The  background  was  assigned  a random  amplitude  and  phase 
value  where  the  amplitude  was  taken  in  the  range  0.05  to  0.2  and  the 
range  of  phases  taken  was  0 to  ir.  Figure  3(a)  illustrates  the  image 
resulting  from  simulating  a receiving  array  with  a beamwidth  of  1.3 
degrees  at  -3  dB  and  a 1-metre  pulse  length.  Figures  3(b)  and  (c)  show 
the  same  image  assuming  increases  by  a factor  of  2 and  4,  respectively, 
in  both  beamwidth  and  pulse  length.  The  target  area  shown  represents  a 
sector  coverage  of  60  degrees.  The  target  is  stretched  in  azimuth  due 
to  the  reduced  effective  aperture  and  broadened  beamwidth  for  off-axis 
beams;  i.e.,  BW  = A / sin  9).  Note  that  interference  effects  in  this 

(Text  Continued  on  Page  8) 
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coherent  sonar  system  can  be  seen  in  Figure  3(a)  even  where  the  target 
is  approximately  eight  "beamwidths"  wide.  Recognition  of  the  target 
as  a square  or  rectangle  is  not  possible  in  Figures  3(b)  and  3(c)  be- 
cause of  the  loss  of  sharp  edge  information,  although  target  detection 
is  obviously  possible. 

The  receiving  and  projector  directional  responses  shown  in  Figure  4 
were  used  in  simulating  the  convolution  model.  Target  information  is 
given  in  the  following  values. 


FIGURE  4.  DIRECTIONAL  RESPONSES 
USED  IN  CONVOLUTION  MODEL 
SIMULATION 


nt  = 64,  m^  = 256 

h = 10  m,  R . = 106.5,  R = 110  m 

min  max 

0 < T < 1.0,  T = 0; 

A n 

background 

n.  P 

parameters 

0.7  < TA  < 1.0,  T = 0; 
A n 

target 

p 

parameters 

The  dimensions  of  each 

of  the 

three  targets  were  30  by  96  cm.  The 
-3  dB  beanrwidth  of  the  receive  direc- 
tional response  was  approximately 
20  cm. 

Figure  5 illustrates  the  effect 
of  using  the  directional  response 
(curve  b)  of  Figure  4 with  range  reso- 
lution of  6,  12,  and  24  cm, 
respectively. 

Figure  6 illustrates  the  effect  of 
using  the  directional  response 
(curve  a)  of  Figure  4 with  range  reso- 
lution of  6 and  12  cm  respectively. 


CONCLUSIONS 


The  effects  of  beamshape  and  pulse  length  can  be  observed  in  the 
target  images.  Subtle  degrading  effects,  difficult  to  express  mathe- 
matically, are  observed  as  range  resolution  is  degraded  and  as  sidelobe 
levels  are  increased.  It  is  conceded  that  the  image  is  "fuzzier"  or 
defocused  as  these  degrading  effects  are  applied.  This  obviously  is 
due  to  the  absence  of  sharp  gradients  in  the  degraded  image.  Though 
these  effects  can  be  observed  subjectively,  it  remains  to  establish 
criteria  for  "adequacy"  of  resolution  for  detection  of  specific  ob- 
ject shapes  in  the  presence  of  the  reverberation. 
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b.  Range  resolution  of  12  cm 


c.  Range  resolution  of  24  cm 


FIGURE  5.  C0NV0LUT1  ODEL  SIMULATION  DIRECTIONAL 
RESPONSE  OF  CURVE  (b)  OF  FIGURE  4 
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Range  resolution  of  6 cm 


b.  Range  resolution  of  12  cm 


FIGURE  6.  CONVOLUTION  MODEL  SIMULATION  DIRECTIONAL 
RESPONSE  OF  CURVE  (a)  OF  FIGURE  4 


